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Suspensions of small anisotropic particles, termed 'rheoscopic fluids', are used for flow visuali- 
sation. By illuminating the fluid with light of three different colours, it is possible to determine 
Poincare indices for vector fields formed by the longest axis of the particles. Because this vector 
field is non-oriented, half-integer Poincare indices are possible, and are observed experimentally. An 
exact solution for the direction vector appears to preclude the existence of topological singularities. 
However, we show that upon averaging over the random initial orientations of particles, singularities 
with half-integer Poincare index appear. We describe their normal forms. 

PACS numbers: 47.57.E-, 47.54.-r, 47.80.Jk 



1. Introduction. Fig.[T^ is a photograph of the surface 
of a randomly stirred fluid. The fluid is a suspension 
of elongated microscopic reflective particles, illuminated 
by red, green and blue lights (R,G,B) coming from three 
different directions, as illustrated in Fig. [TJ). Two circles 
indicate points with the property that, on traversing a 
small clockwise circuit around this point, we encounter 
the primary colours without repetition (in one case R- 
G-B, in the other R-B-G). At first sight the existence 
of such points seems unremarkable, but we argue below 
that it is indicative of a singularity in the direction field 
of the particles which appears to be forbidden. Our letter 
explains why such singularities are nevertheless observed, 
and classifies their normal forms. 

The fluid used to produce Fig. [1^ was a commercially 
available rheoscopic fluid, which is used (diluted in water) 
for flow visualisation [1]. It contains microscopic rod-like 
crystals. It is desirable to understand what aspects of the 
fluid motion are revealed by the rheoscopic fluid (various 
aspects are discussed in [H, 0, S 0] ) • In order to pose 
the questions in the simplest setting, we consider incom- 
pressible flow in two dimensions, before finally consider- 
ing the three-dimensional case. The reflective elements 
are assumed to be rod-like axisymmetric crystals which 
are very short compared to the characteristic length scale 
of the velocity field of the flow, v{r^t). 

In the following we show that Fig. [1^ demonstrates 
the existence of singularities in the direction field which 
have a half-integer Poincare index. (The Poincare index 
is a topological invariant. For a vector field in the plane, 
the Poincare index of a closed curve is the number of 27r 
clockwise rotations of the vector field as the curve is tra- 
versed, also clockwise. Curves with a non-zero Poincare 
index encircle a singularity of the field.) 

There are two reasons why the appearance of such sin- 
gularities is unexpected. Intuitively, it might be expected 
that the long axis of the particles will align with stream- 




FIG. 1: (a) Shows textures obtained by illuminating the sur- 
face of randomly stirred rheoscopic fluid with diffuse red, 
green, and blue light sources, as illustrated in (b). The circles 
in (a) indicate the positions of singularities. 



lines of the flow. Streamlines around a vortex or about a 
saddle point have Poincare index of ±1. However, singu- 
larities of this type would result in the particles reflecting 
each colour twice upon making a circuit (see Fig. [2^), 
which is not consistent with Fig. [1^. A more persua- 
sive argument is that the exact solution (given below) of 
the equation of motion for the axis of the crystals shows 
that the Poincare index of this vector field is invariant, 
which seems to preclude the emergence of patterns such 
that shown in Fig. [1^. However, to compute the inten- 
sity of reflected light one must average over the random 
initial particle orientations. We show that this gives rise 
to an order-parameter field which does exhibit singulari- 
ties. We remark that, while the use of colours to enhance 
rheoscopic visualisations was previously demonstrated in 
[7], the technique was not used to reveal singularities. 

The axial direction field of the rod-like particles is 
non-orientable (that is, the sign of the vector is irrel- 
evant). This allows other types of singularity, such as 
those shown in Figs. [2)3 and c, which have a half-integer 
Poincare index, and which are consistent with the par- 
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FIG. 2: (a) Shows the direction field around a vortex with 
Poincare index 1 and the corresponding colour texture cycling 
through R ^ G ^ B ^ R ^ G ^ B. (b), (c) Show the same 
for a core singularity with index ^ and for a delta singularity 
with index — | respectively. 



tides reflecting each colour only once upon traversing 
a closed curve. These singularities have not previously 
been considered in fluid dynamical problems, although 
they are seen in fingerprints (where they are known as 
the core and delta, respectively [8]). Fig. [1^1 is evidence 
that they are present in our experiment. 

2. Equation of motion and its solution. The equa- 
tion of motion for an ellipsoid of revolution in a flow at 
low Reynolds number was obtained by Jeffrey [9], and 
Bretherton \ld\ showed that the same equation applies 
to a general axisymmetric body. In the limit of rod-like 
particles, a unit vector n{t) in the direction of the axis 
of the microscopic particles satisfies 



h — An — (n • An)n 



(1) 



where A(r(t), t) is the velocity-gradient tensor at the cen- 
tre r{t) of the particle, with elements Aij = dvi/drj. The 
particle position advected by the flow, r = v{r,t). 

An exact solution of the equation of motion ([T]) is ob- 
tained from a matrix M(t) which is obtained by integrat- 
ing the linear differential equation 



M = A(r(t),t)M 



(2) 



where r{t) is the trajectory of the centre of the rod. This 
matrix is the monodromy matrix describing the evolu- 
tion of the infinitesimal separation Sr{t) of neighbour- 
ing points in the flow: we have Sr{t) = M(t, to)^r(to), 
where Sr{to) is the initial infinitesimal separation at time 
to. The initial condition for equation ([2j) at time to is 
M(t, to) = I, where I is the identity matrix. Now if no 
is the initial direction of the rod at time to, the direction 
at time t is given by 



n{t) = d{t)/\d{t)\ , d{t) = M(t,to)no . 



(3) 



This solution was first given by Szeri [4]. We are inter- 
ested in the vector field of the rod orientations at position 
r and time t. The matrix M then depends upon position 
as well as time, and we write M(r,t,to) for the mon- 
odromy matrix of a trajectory which reaches r at time t. 



starting from ro at time to- It follows from (|3]) that the 
vector field of rod orientations is 



n(r,t) = 



M(r,t,to)no(ro) 
|M(r,t,to)no(ro)| 



(4) 



The monodromy matrix M(r,t,to) is a smooth func- 
tion of the final position of the trajectory, r. The solution 
(IH) can therefore only be discontinuous if the initial direc- 
tion field is discontinuous, or if the denominator |Mno| is 
equal to zero, which is not possible because det(M) = 1. 
If the initial direction vector field no(r) is non-singular, 
we therefore conclude that the direction field n(r,t) re- 
mains non-singular for all times. Because the vector field 
generated by (|4j) is smooth, the Poincare index of this 
field is zero for any closed curve. 

In an earlier paper [11] we discussed the textures 
formed by rods-like particles in a complex flow for a spec- 
ified initial direction field no(r). Here we deal with the 
more complex case where we must average over the ran- 
dom initial orientation. We show below that this leads 
to singularities of the orientation field. 

3. An order parameter for rheoscopic Euids. Initially, 
at time to, the rod-like particles in a rheoscopic fluid 
are randomly oriented, due to the effects of Brownian 
motion. We must therefore consider the distribution 
of rod directions generated by the solution (|4]) at each 
point in the flow. This can be described by a probability 
density P{6) for the rod orientation angle 6 (satisfying 
P{0 -h tt) = P{0), because the rods are non-oriented). 
This probability density depends upon both position and 
time. Specifying the angle distribution at each point in 
space and time would provide too much information to 
be a useful description. It is therefore desirable to map 
this distribution P(^,r,t) to an order parameter vector 
(^(r,t). The direction of this non-oriented vector should 
represent the predominant direction of the rods, and its 
magnitude should indicate the degree of ordering (with 
|(^| = 1 when the rods are all in the same direction, and 
ICI = when their angular distribution is isotropic). 

The initial direction vector no in ([2j) is a random vec- 
tor, uniformly distributed about the unit circle. The vec- 
tor d{t) = M(t,to)no is therefore distributed about an 
ellipse. If the unit circle is represented as X ' X — 1, the 
ellipse upon which d lies is represented by the equation 



x-Kx = l, K=(MM^)- 



(5) 



This ellipse has its longest axis aligned along a direction 
0, which is in the direction of the eigenvector correspond- 
ing to the largest eigenvalue of M M^. Its aspect ratio u 
is the square root of the ratio of the eigenvalues of MM^ 
(we choose to consider ly > 1). It is natural to define the 
magnitude of the order parameter to be a function of u 
which interpolates between zero (when u = 1) and unity 
(as ^ oo). In a later paper we shall argue that the fol- 
lowing expression is the most natural definition for the 
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order parameter: 

C=^n(.-), (6) 

where n(^) is a non-oriented unit vector with angle 9. 

4. Relation between order parameter and light scatter- 
ing. The colour which is reflected by the rheoscopic fluid 
in the experiment illustrated in Fig. [1] may be related 
to the order parameter. The details of this relation de- 
pend on the ratio of the length of the rod-like particles in 
the rheoscopic fluid to the wavelength of light, and upon 
their surface roughness. For illustration we discuss the 
simplest case, where the rods are short compared to the 
wavelength of light. In this limit the scattered intensity 
from a rod is proportional to the square of the projected 
cross-section, so that a rod at angle 6 scatters light from 
a source which is perpendicular to the direction (j) with 
an intensity proportional to cos^(^ — (f)). The combined 
effect of scattering from the three light sources results in 
additive colour mixing, so that at any given point the re- 
flected colour C is a weighted combination of red, green 
and blue (i?, G, B) of the form 

C = /(O) R + /(27r/3) G + /(47r/3) B (7) 

where is the average of cos^(6> — (j)) over the rod 

orientation: 

n27T 

I {(I)) = dO P{0) cos2(0 - 0) . (8) 
Jo 

After a lengthy but elementary calculation we obtained 
the following expression for the reflected colour in terms 
of the parameters ((9, v) of the order parameter: 

C = — ^ — -\4R[iycos^0^sin^0] (9) 

+ G[2{l-u) cos^ 6 - 2\fz{y- 1) sin OcosO^ 3i/+ 1] 
+ 5[2(l-i/)cos^^+2y3(i/-l)sin^cos^ + 3i^ + l]} . 

Figs. [3^1 and c show the order parameter obtained from 
dH) and (|6]) for small rods moving in a random flow field 
(we used the same model as in [uj), starting from a uni- 
form distribution of angles. In the vicinity of zeros of the 
order parameter, there are patterns which resemble the 
delta (Fig. [3ti) and core (Fig. [3t) singularities illustrated 
in Figs. [2t and b. 

The remaining panels of Fig.[3](b and d) show the light 
reflected from three coloured sources, computed using 
equation (|9|). (The additive colour mixing for this figure 
was performed using MatLab.). These singularities have 
a Poincare index of and they are therefore consistent 
with the experimental result shown in Fig. [1^. In the 
remainder of this letter we discuss how the occurrence of 
these structures can be understood, and the possibility 
that additional types of singularity might appear if the 
flow is three-dimensional. 




FIG. 3: (a) Shows the order parameter field for a random 
flow, (b) Shows the same image colour coded using equation 
(jS)) to indicate the additive colour mixing of light from three 
different sources. The circle marks the position of a delta 
singularity; (c) and (d) the same, but for a core singularity. 

5. Singularities of the order parameter. If zeros of the 
vector order parameter field exist, the Poincare index of a 
curve may be non-zero even though the order parameter 
depends smoothly upon position. A singularity where the 
order parameter is equal to zero occurs where the ellipse 
upon which the vector d{t) = M(t, to) no lies degenerates 
to a circle. Thus zeros of the order parameter occur when 
M is a rotation matrix. 

First we consider whether such singularities are 
generic. Because the two-dimensional flow is area- 
preserving, the 2x2 matrix M satisfies det(M) = 1. 
This matrix can be written in a form determined by three 
parameters A, hz^ x- 

M = D(A,A-i)S(/t)0(x) (10) 

where D, S and O are respectively diagonal, shear and 
rotation matrices: 

D(A,.A,, = «) , SM=(J^), 

y— smx cosx J 

The singularity occurs when A = 1 and = in (p!Q|) 
(with no condition upon x), which is realised upon vary- 
ing two parameters. Since the coordinate space is two- 
dimensional, the singularities occur at isolated points in 
the plane. 



Having identified the condition defining the singular 
point, we now turn to consider the form of the order pa- 
rameter field in its vicinity. To leading order, generically 
the parameters A and hc depend linearly on position in the 
vicinity of a singular point at which the monodromy ma- 
trix becomes a pure rotation. Let us assume that there 
is a singular point at vq = (xq, yo). In the vicinity of this 
point there exists a coordinate system R = (X, Y) such 
that the monodromy matrix is in the normal form 

M{R) = D(l + iX, S(±r) 0{x)+OiR^) . (12) 

The local coordinate system is related to r by a linear 
transformation: r — tq = Ti^, where T is a 2 x 2 matrix 
and the sign in (p!2|) is chosen so that T is a non-inverting 
transformation (det(T) > 0). The vector order param- 
eter for the normal form is plotted in Fig. [H for both 
choices of the sign in (p!2|). and it can be seen that these 
two patterns resemble the core and delta singularities of 
fingerprints (as described in Q). It is instructive to anal- 
yse the behaviour of the order parameter in the vicinity 
of these singularities. Retaining only the leading order 
terms in X and Y from ([T2l). we find 



1 + X ±Y 
±Y 1-X 



■0{R^) (13) 



where R = \J + Y^ . The eigenvalues of MM^ are 
A± = l±i? + 0(i?2), so thatz/= v^A+ZA- = R^O{R^). 
Writing (X, F) = (i? cos 6, sin 6) and (cos <I>, sin <l>)^ 
for the eigenvector of MM^ corresponding to A+, we 
find <I> = ±^0, so that 



C{R,e) = Rn{±le) ^ 0{R^) 



(14) 



If the positive sign is chosen, points in the radial di- 
rection (<l> = Gmodyr) for only one ray (0 = 0). For the 
negative sign, points radially outwards along three rays 
(O = 0, ±27r/3). These properties are characteristic of 
the core and delta singularities, respectively. 

6. Three-dimensional Eows. Thus far, we have consid- 
ered two-dimensional flows. In the experiment illustrated 
in Fig.[T^, the depth of the rheoscopic fluid was a few mil- 
limetres, and comparable to the scale size of structures in 
the texture shown in Fig.[T^, which shows a region which 
is a few centimetres across. We should therefore consider 
what additional structures might arise when the flow is 
three-dimensional. In a three dimensional flow the di- 
rection vector n(r, t) covers a sphere rather than a circle 
(it remains non-oriented). The optical depth was very 
small, so that light was reflected from a thin layer of 
fluid just below the surface. It is only the direction of 
the projection of n(r, t) in the plane of the fluid surface 
which determines the colour of the scattered light. This 
projected vector field, np(r,t), can have simple zeros if 
there are positions where the rods point out of the surface 
of the liquid. The projected vector field np(r, t) then has 




FIG. 4: (a) Shows the order parameter field C(-^) for the 
canonical singularity of eq. (|12|) . with a positive sign multi- 
plying Y. This is a core singularity, (b) Shows the order 
parameter field when the negative sign is chosen. This is a 
delta singularity. 



a simple zero, with Poincare index +1. We conclude that 
if the flow is three dimensional, double rotations of the 
primary colours (illustrated in Fig.[2ti) may be observed, 
as well as single rotations (as in Figs. [2)3 and c). 

7. Concluding remarks. We have shown that our sim- 
ple experiment on a rheoscopic fluid reveals singularities 
of the direction field with Poincare index equal to 
We defined an order parameter for this field, and showed 
that its zeros have normal forms which are analogous to 
the core and delta singularities of fingerprints. The dis- 
tinctive behaviour of the order parameter in steady flows 
will be addressed in a longer paper, in preparation. 
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